www.gusucode.com > 小波分析理论与Matlab 7实现 > 小波分析理论与Matlab 7实现/code/Ch16/examp16_23.m

    %装载并显示原始图像
load bust;
%加入含噪
init=2055615866;
randn('seed',init);
X1=X+20*randn(size(X));
subplot(221);image(X1);
%colormap(map);
title('原始图像');axis square;
%用小波db4对图像X进行一层小波包分解
T=wpdec2(X1,1,'db4');
%重构图像近似部分
A=wprcoef(T,[1 0]);
subplot(222);image(A);
title('图像的近似部分');
axis square;
%边缘检测
%%原图像的边缘检测
BW1 = edge(A,'sobel');
subplot(223);imshow(BW1);
title('原图像的边缘');
axis square;
%图像近似部分的边缘检测
BW2= edge(X1,'sobel');
subplot(224);imshow(BW2);
title('图像近似部分的边缘');
axis square;